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We show that perturbation theory may give reasonable numbers for the 
decays of the bottomonium and charmonium ground states to e~'"e~ and to 
77. To reach this conclusion it is important to perform the resummation of 
logs. In particular, we obtain the value r(77b(15) ^77) = 0.35 ± 0.1(th.) ± 
0.05(Aqcd) KeV. 

1. Introduction 

In Ref. yy, the decays of Heavy Quarkonium to e+e^ and 77 were 
computed with next to leading log (NLL) accuracy within perturbation 
theory. When presenting these results in this conference, one (reasonable) 
complaint was the complete absence of numbers. We would like to fill this 
gap by doing the phenomenological analysis of these results that we quote 
here for ease of reference. 

We first quote the matching coefficients at the hard scale |3 El • 
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where Cf = {N^ — l)/{2Nc), whereas the renormalization group improved 
matching coefficients at NLL (for the vector and pseudo-scalar) read (m^) 
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where /3o = ^-Ca — ^TpUf, z = 
coefficients Ai in Eq. Q read 

7rC/[3/?o(26C72 + 19C7aC/ - 32Cj) - Ca(208C^ + 651CaC/ + IIGC^)] 



^3 

A, 



78 (3iC A 

PoiMs + 1) - 3) + Ca(15 - Us{s + 1)) 



6(/3o - 2C^)2 
247rC2(3/3o - 11Ca){5Ca + 8C/) 



13Ca(6/3o-13C^)2 



/32(6/3o-13C7a)(/3o-2Ca) 



C72 (-9Ca + lOOC/) 
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+PoCA{-7ACf + C7a(42 - 13s(s + 1))) + 6/3o(2C/ + Ca(-3 + s(s + 1))) 

By setting Vp ~ mCfas/n, Bs{vp) includes all the large logs at NLL order 
in any (inclusive enough) S-wave heavy-quarkonium production observable 
we can think of. For instance, the decays to e^e~ and to two photons at 
NLL order read 



T{VQ{nS)^e+e- 



r(PQ(n5)^77) = 2C7a 
^ 2Ca 




n 



(5) 



mqCfa^ 



n 



mqCfas 



n 



n 2 



PQ(nS)_ 



n 



^\l + 2{Bi{vp)-l) + 25^n} , 
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{l + 2(SoK)-l) + 2<50„} 



where V and P stand for the vector and pseudo-scalar heavy quarkonium, 
as = tts(^'p)) and (^'^(z) = ^ ^^Jn^^^ and T(z) is the Euler F-function) 
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which has been read from Ref. 

It is not our aim to perform a full fledge analysis of Eqs. © and © here 
but rather to see what are the general trends obtained by the introduction 
of the resummation of logs, as well as to give some predictions when the 
results appear to be reliable enough. 

2. Phenomenological analysis 

In this section we perform the phenomenological analysis of the above 
results for the bottomonium and charmonium systems. 

2.1. b-b IS states 

We first consider the b-h IS states and their decays r(T(15) — > e+e~) 
and r(%(15) 77). For the first decay we will be able to compare our re- 
sults with experiment whereas our numbers for the second can be considered 
to be a prediction. 

We plot r(T(lS') e^e~) in Fig. ^ versus the renormalization scale 
u. We consider the LO/LL result (they are equal), the NLO result and the 
NLL result. We can see that the LO/LL result can match the experimental 
figure for a reasonable value of i^, close to the soft scale. This agreement 
is destroyed once the NLO is considered. The best value we can obtain is 
0.667 KeV. The reason seems to be the fact that we now have two scales in 
the problem: the hard (~ m) and the soft (~ mas) scale. The final outcome 
is that perturbation theory breaks down before the normalization scale v 
can reach the typical soft scale of the problem. The standard solution to 
this problem goes by doing a renormalization group analysis, summing up 
all the large logs that appear in the problem. Actually, in our case, we are 
going to have large logs produced by the ratio of the hard and soft scale and 
by the ratio of the soft and ultrasoft scale (~ ma^). We can see that the 
use of the NLL expressions improves the agreement with the experimental 
result (the best value now becomes 0.837 KeV) and enlarges the range of 
applicability of perturbation theory. Moreover, the expansion seems to be 
convergent, being the NLL result a correction with respect the LL order 
one. Nevertheless, perturbation theory still seems to break down before the 
renormalization scale u can reach the typical soft scale (although getting 
closer to it than in the NLO calculation) and the optimal result is off the 
experimental value by around 50%. This is a large effect. Therefore, in 
order to confirm this picture, a full NNLL result should be obtained. This 
is a difficult computation but were the convergent pattern confirmed the 
outcomes will be important. One could then start quantifying the size of 
the non-perturbative effects in a reliable manner. 
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Fig. 1. Plot o/r(T(lS') e+e") with LO/LL (dashed line), NLO (dot-dashed line) 
and NLL (dotted line) accuracy versus the renormalization scale v. The horizontal 
line and its hand give the experimental value and its errors: r(T(lS') e+e^) — 
1.314 ±0.029 KeV\^. 



We now perform a similar analysis for r(%(15) — > 77). For this decay 
no experimental figure exists. Therefore, we will be able to give a predic- 
tion for it. The picture is completely analogous to the r(T(15) e~^e~) 
case. Actually the scale of minimal sensitive is a little bit smaller than in 
the previous case, which makes us more confident about the result. This 
confidence will be further boost below when we consider the charmonium 
case for which experimental numbers exist. We anticipate that very good 
agreement is reached in that case. We then give a prediction for this decay: 

T{rib{lS) ^ 77) = 0.35 ± 0.1(th.) ± 0.05(AQCD)KeV, (8) 

where "Aqcd" stands for the error due to the variation of Os {as{Mz) = 
0.118 lb 0.003) and "th" for the theoretical errors. The later are difficult to 
obtain and here we only pretend to roughly estimate their size. We con- 
sider the variation of the decay if we increase v by two GeV with respect 
the optimal scale. This gives the theoretical error quoted in Eq. ©. We 
consider this estimate conservative in view of the good agreement with data 
obtained for the r/c case. Note that to have larger errors would be incon- 
sistent with the assumption that the remaining corrections (perturbative 
and non-perturbative) are smaller, or at least of the same order, than the 
difference between the LL and NLL result. 
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Fig. 2. Plot ofT[rib{lS) -> 77) with LO/LL (dashed line), NLO (dot-dashed line) 
and NLL (dotted line) accuracy versus the renormalization scale v. 



For the above states one may consider reasonable to believe that a per- 
turbative approach is a good starting point for these systems, since the 
soft scale is clearly in the perturbative regime. Nevertheless, one should be 
careful since the ultrasoft scale also enters the game through the matching 
coefficient (ultrasoft effects first appear at 0(af In^Os)). This means that 
for V < 1.5 GeV, as{iy'^/mi,) starts to blow up. Obviously, the situation 
is even worse for the n = 2 bottomonium states, since the soft scale goes 
divided by 1/n? (actually, we anticipate that a bad behavior is found in 
this case). Surprisingly, however, for the charmonium system we are in a 
situation quite similar to the n = 1 bottomonium states since, although the 
typical soft scale is smaller than the n = 1 bottomonium soft scale, this is 
compensated by the fact that the charm mass is smaller than the bottom 
one, so that we can run down u to scales quite close to the typical soft scales 
of the problem before the break down of perturbation theory takes place. 
With all these qualifications let us consider the n = 1 charmonium or n = 2 
bottomonium states and see whether reasonable numbers are obtained. 



2.2. n = 1 charmonium states 

We study r(J/^'(15') e~^e~) in Fig. ^ Surprisingly a pretty similar 
picture than in the r(T(lS') — > e~^e~) appears. Actually, the difference with 
experiment is of the same order, around 50 %. 

We now consider T{rjc{lS) 77) in Fig. 0] In this case we have exper- 
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Fig. 3. Plot of r{J/<i/{lS) e+e") with LO/LL (dashed line), NLO (dot- 
dashed line) and NLL (dotted line) accuracy versus the renormalization scale 
V. The horizontal line and its hand give the experimental value and its errors: 
T{J/^{\S) e+e-) = 5.14 ± 0.31 KeV J/. 



imental data to compare with, unlike the r/^ case. A similar pattern to the 
r/6 case is found and, moreover, we get agreement with experiment. This 
is quite shocking since, at the scale of minimal sensitive, the perturbative 
result hits the experimental one just in the middle. Note that to get this 
agreement the resummation of logs is necessary. The scale of minimal sensi- 
tive is around 850 MeV, of the same order than the typical soft scale of the 
state. However, these numbers should be taken with caution at this stage, 
since the running involves the ultrasoft scale, which has been run down to 
very low scales ~ 500 MeV in a correlated way. 



2.3. n = 2 bottomonium states 

We now turn to the n = 2 bottomonium states. We show our results in 
Figs. Elandini In principle, this is the most problematic case since the soft 
scale of the T{2S) is ~ l/n^(= 1/4) x the soft scale of the T(IS') (even 
if partially corrected by the fact that Os would be larger for the 2S state 
than for the IS state). Actually for r(T(2S') e^e~) we can compare with 
experiment and our result is a factor two smaller than the experimental 
number. We note that, even for the leading order result, we have to go to 
very small scales to get agreement with experiment (~ 600 MeV). This raises 
doubts about our perturbative analysis for the n = 2 bottomonium states. 
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Fig. 4. Plot of r{r]c{lS) 77) with LO/LL (dashed line), NLO (dot-dashed line) 
and NLL (solid line) accuracy versus the renormalization scale v. The horizontal 
line and its band give the experimental value and its errors: T{r]c{lS) 77) — 
7.2 ± 1.2 KeVlBl- 



Somewhat, even if the resummation of logs helps, perturbation theory still 
breaks down before one can reach the typical soft scale of the problem (which 
is very small). 



2.4- Final discussion 

Even if one should wait for the complete NNLL result, we can start 
to see some trends. The uncalculated (perturbative or non-perturbative) 
terms seems to be larger for the decays into e~^e~ than for the decays into 
two photons. One could then start to speculate about the size of the non- 
perturbative effects in each case. Good enough, the agreement obtained 
with the experimental figure of T(r]c{lS) — > 77) makes us quite confident 
that the result we have obtained for T{7]h{lS) 77) will be quite close to 
the experimental number. 



3. Conclusions 

We have performed a phenomenological analysis of the NLL results ob- 
tained in Ref. for the heavy quarkonium decays to e'^e~ and to two 
photons. 
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Fig. 5. Plot ofr{T{2S) e+e") with LO/LL (dashed line), NLO (dot-dashed line) 
and NLL (dotted line) accuracy versus the renormalization scale v. The horizontal 
line and its hand give the experimental value and its errors: r(T(25) e+e^) — 
0.576 ±0.024 KeV\^. 
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Fig. 6. Plot of T{riii(2S) 77) with LO/LL (dashed line), NLO (dot-dashed line) 
and NLL (dotted line) accuracy versus the renormalization scale v. 



For r(T(15) e^e ), we find that the resummation of logs signifi- 
cantly improves the agreement of the perturbative result with experiment 
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in relation with a pure NLO evaluation, such that, at the place of minimal 
sensitive, the theoretical number is off the experimental one by around 50%. 
Moreover, overall convergence is found, with the NLL order result being a 
correction with respect the LL order one. Note that in order to be so, we 
have to take v of the order of the soft scale such that the large logs are re- 
summed. Surprisingly, the very same picture holds for T[J /"^{IS) e+e~) 
too. It would certainly be challenging to have the complete NNLL result to 
check whether this pattern of convergence survives, since the difference with 
experiment is still large in both cases. If this is so one could start to reliable 
estimate non-perturbative effects in heavy quarkonium (and its relation with 
the ultrasoft cutoff) and, for the charmonium case, support the view that 
one can actually use perturbation theory as an starting point for its study 
(see 0). On the other hand, for r(T(25') — > e'*"e~), we find an strong dis- 
crepancy with the experimental figure, rasing doubts that one can actually 
use perturbation theory there. This could make sense since, for Coulomb- 
type bound states, the soft scale of the T{2S) would be ~ l/n?{= 1/4) x 
the soft scale of the T{\S) (partially corrected by the fact that would 
be larger for the 2S state than for the IS one). Even if the resummation of 
logs helps, perturbation theory still breaks down before the renormalization 
scale v can reach the typical soft scale of the problem. Somewhat we feel 
that something similar still happens for the T(15') and the J /"^{IS) in a 
less severe way, since in these last cases one can get much closer to what we 
believe are the typical soft scales of the problem. In any case, we would not 
like to draw definite conclusions just from this analysis. 

For the decay of the heavy quarkonium to two photons we only have 
experimental data for the Charmonium case. In this case perfect agreement 
with experiment is obtained. This boosts our confidence that our number 
for r(r/fe(15) — s- 77) will be reliable, which we take as one of the major 
results of this paper. 

For T{r}i,{2S) — > 77), we prefer not to draw any conclusion in view of 
the failure of the results obtained for r(T(25) e^e~). 

It is a general trend that the resummation of logs improves the agreement 
with the experimental result (when available). This can be understood by 
the fact that makes perturbation theory stable up to smaller scales. 

We can also see that the picture is quite similar to the one obtained in 
Ref. [SI for the hyperfine of the heavy quarkonium, which has recently been 
computed with NLL accuracy. 

The fact that we obtain reasonable numbers for the charmonium sys- 
tem can be considered a surprise. One may think that in this case one has 
run down to very low scales the ultrasoft scale. Actually, at the numerical 
level, the ultrasoft scales for the ground states of the bottomonium and 
charmonium seem to be similar. The reason is that, even if the soft scale 



10 



USTRON PRINTED ON FEBRUARY 1, 2008 



of charmonium is smaller than of bottomonium, this is compensated by the 
fact that the charm mass is smaller than the bottom mass. This explains 
the similar behavior found for both systems with errors of the same size. 
For the n = 2 bottomonium states, the behavior seems to be different, with 
a major breakdown of perturbation theory, making the results not trustwor- 
thy (actually a factor two disagreement with experiment is obtained when 
experimental results are available). In any case, even for the bottomonium 
and charmonium ground state, the ultrasoft scale has been run down to 
very low scales. This is a potential problem of the whole analysis. The issue 
would be certainly clarified if a complete NNLL computation were available 
for the decays. In case a convergent pattern is observed for the perturba- 
tive series, it would certainly be considered an indication that perturbation 
theory can be applied for these systems. One can then start considering 
a quantitative study of non-perturbative effects and its relation with the 
ultrasoft cutoff. In particular one may start to consider what the ratio be- 
tween Aqcd and mv'^ exactly is and try to apply the results obtained in 
Ref. PI for the non-perturbative corrections. 
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